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In this letter we study the Haldane model with on-site repulsive interactions at half-filling. We 
show that the mean-field Hamiltonian with magnetic order effectively modifies parameters in the 
Haldane Hamiltonian, such as sublattice energy difference and phase in next nearest hopping. As 
interaction increases, increasing of magnetic order corresponds to varying these parameters and con¬ 
sequently, drives topological transitions. At the mean-field level, one scenario is that the magnetic 
order continuously increases, and inevitably, the fermion gap closes at the topological transition 
point with nonzero magnetic order. Beyond the mean-field, interaction between fermions mediated 
by spin-wave fluctuations can further open up the gap, rendering a first-order transition. Another 
scenario is a first-order transition at mean-field level across which a canted magnetic order develops 
discontinuously, avoiding the fermion gap closing. We find that both scenarios exist in the phase dia¬ 
gram of the Haldane-Hnbbard model. Our predication is relevant to recent experimental realization 
of the Haldane model in cold atom system. 


Correlation and topology are two of central topics in 
modern condensed matter physics. The Hubbard model 
of spin-1/2 fermions with on-site interaction is one of 
the most famous model that gives rise to highly non¬ 
trivial correlation effects, such as antiferromangetic order 
at half-filling. This model is relevant to many strongly 
correlated materials such as High-Tj, cuprate supercon¬ 
ductors. Recently, this model has also been simulated by 
using ultracold fermions in optical lattices [IHS]. In 1988, 
Haldane proposed a model of noninteracting fermions in 
the honeycomb lattice, which can give rise to topological 
band structure and quantized Hall conductance without 
external magnetic field [3]. This effect is now known as 
the quantum anomalous Hall effect, and has been ob¬ 
served experimentally in magnetically doped topological 
insulators [S] . The Haldane model has also been realized 
recently in cold atom experiment using shaking lattice 
technique 0111. 

Motivated by recent cold atom realization of the Hal¬ 
dane model, in this letter we study the Haldane-Hubbard 
(HH) model of spin-1/2 fermions. In this model each spin 
component experiences the same single-particle Hamil¬ 
tonian described by the Haldane model. At half-filling 
(A/ = = Ns/2, Ns is the number of sites), the sys¬ 

tem is a topological band insulator in the noninteracting 
limit. We consider only on-site repulsive interaction be¬ 
tween two spin components as in the Hubbard model. 
In the strongly interacting limit, the system will enter a 
Mott insulator phase. Therefore, as interaction increases, 
we expect that two things will happen: One is the devel¬ 
opment of certain magnetic order, and the other is the 
transition from a topological band insulator to a topolog¬ 
ically trivial Mott insulator. A natural question is how 
these two phenomena influence each other. The study 
of this question will shed light on the interplay between 
correlation and topology [8j. Previously, although there 
have been considerable interests focusing on the Kane- 
Mele-Hubbard model 0, the HH model is much less in¬ 


vestigated m- 

In this letter we show that, at the mean-field level the 
magnetic order can drive a topological transition, either 
by a continuous second-order transition with fermion gap 
closed at finite magnetic order, or by a first-order tran¬ 
sition with a jump of magnetic order parameter. Within 
the mean-field calculations, we show that both scenar¬ 
ios exist in the phase diagram of the HH model. We 
then go beyond the mean field by considering fluctua¬ 
tions of magnetic order. We show that for the former 
case, the spin-wave fluctuations generate effective inter¬ 
action among gapless fermions at the nominal critical 
point of topological transition. This interaction, if suffi¬ 
ciently strong, can open up a gap and drives the transi¬ 
tion to first-order. 

The Model. We consider the HH model on a honey¬ 
comb lattice whose Hamiltonian is given by 

.Hhh = Hii + U ^ (1) 

i 

[clsCj,s + h.c.) 

(b>:S 

-t2 X (2) 

{{ij)),s i,s 

where s = ± refers to spin up and down, respectively, 
the t 2 “term represents next nearest hopping with a non¬ 
trivial phase (/ij = ±(/ for different sublattices, and the 
M-term adds a potential imbalance between A and B 
sub lattices, as = ±1 for i belonging to A or H sublat¬ 
tices. The t 2 term breaks time-reversal symmetry, and 
the Af-term breaks the inversion symmetry. Both terms 
open up the gap at Dirac points, and for half-filling, a 
phase diagram (without the interaction term) including 
a topological transition from trivial insulator to topolog¬ 
ically nontrivial insulator is shown in Fig. 1(a), across 
which the gap at one of the Dirac point is closed. For 
the topologically nontrivial insulator, each spin compo- 
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FIG. 1: (a) Phase diagram for the noninteracting Haldane 
Hamiltonian Hyi (Eq. 1^, in terms of AB sublattice energy 
imbalance M and phase (f). The three trajectories labelled by 
a, b and c correspond to the evolution of mean-field Hamil¬ 
tonian Hmf (Eq. 1^ with the increasing of magnetic order, 
as U/t\ increases, for different 12 /ti as labelled by (a-c) in 
Fig. H The solid line means continuous evolution while the 
dashed line means discontinuous jump. Panels (b) and (c) 
illustrate two types of AF magnetic order: (b) is collinear AF 
order and (c) is canted AF order with nonzero spin chirality. 

nent fills the lower band with Chern number equalling 1 
and the total Chern number C = 2. The interaction term 
can be decoupled as 

i i 

+ + (3) 

i ^ ' 

where we have introduced an on-site magnetic order pa¬ 
rameter 417(S.i)/3 = m^. Thus the mean-field Hamilto¬ 
nian Hmf is given by 

Hmf = Hfl — ^2 

i 

Relation between AF order and Topology. Before pro¬ 
ceeding to the self-consistent mean-field calculation, we 
would like to first discuss the relation between the follow¬ 
ing two types of possible AF order and the parameters 
in the Haldane model. 

(A) Collinear AF order, i.e. nii = m on the sublat¬ 
tice A and m.i = — m on the sublattice B, as shown in 
Fig. [^b). Due to the spin rotational symmetry, we can 
always choose m = mz. Thus, it adds a spin-dependent 
contribution on M in the single-particle Haldane model 
of Eq. i.e. M ^ M -\- sm, where s = ± denotes spin. 



FIG. 2: Mean-field phase diagram for HH model in the 
{t 2 jti,U/ti) plane. The green solid line is a second-order 
transition between a topological band insulator (Phase I) 
and a topological collinear AF insulator (Phase H). The blue 
dashed line is a topological transition between Phase H and a 
topological trivial collinear AF insulator (Phase HI). The red 
dash-dot line is a first-order transition between either Phase 
I or Phase H and a topological trivial canted AF insulator 
(Phase IV). The blue solid line is a second-order phase tran¬ 
sition between Phase HI and Phase IV. Here (a-c) label three 
trajectories discussed in Fig. [^and Fig. 


(B) Canted AF order. For simplicity, we consider the 
situation that m are different among the three A sites 
(denoted by i, j and k) of one honeycomb, as shown in 
Fig. 1(c), which leads to finite “scalar spin chirality” 
order S = (Si) ■ ((Sj) x (S^)) ( Since we are concerned 
with magnetically ordered state, we do not use the usual 
definition, i.e. (Si • (Sj x S^)). ). Within each unit cell, 
tiii at A site is approximately opposite to mi at B site. 

Then we can apply an on-site spin rotation Ui so 
that Uj{mi ■ Si)Ui = ^Imijsi^ for A sublattices and 
Uj{ini-Si)Ui = — ||mi|si 2 , where Siz is the Pauli matrix 
associated to spin. Qualitatively speaking, this local spin 
rotation introduces an additional Berry phase factor ±(j) 
in the next nearest hopping term for different sublattices, 
where cj) is approximately one sixth of the solid angle ex¬ 
panded by mi,mj,mfc, that is to say, (j) in the original 
Haldane Hamiltonian Eq. should be replaced by an 
effective phase 4>eS = (f> + 4>. 

Therefore, the mean-field Hamiltonian with a collinear 
AF order corresponds to a free-fermion Hamiltonian Hu 
with a modified spin-dependent effective M, and the 
mean-field Hamiltonian with a canted AF order corre¬ 
sponds to a free-fermion Hamiltonian with both M and (j) 
in i?H modified. Thus, the noninteracting phase diagram 
in Fig. I^a) is helpful for understanding the mean-field 
phases, with M and (j) replaced by effective parameters 
determined by magnetic order. There emerge two differ¬ 
ent scenarios about how magnetic orders drive transition 
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between topological nontrivial and trivial insulators. 

(i) If AF order increases continuously as interaction 
strength U increases, the mean-field Hamiltonian will 
evolve continuously cross the phase boundary from C = 2 
to C = 0 insulator. Inevitably, there will be a topologi¬ 
cal transition at which gapless fermions and a finite AF 
order coexist. 

(ii) A first-order transition occurs as U increases, at 
which a jump of AF order brings the system from C = 2 
regime in the phase diagram to C = 0 regime. 

Mean-field Phase Diagram. A self-consistent mean- 
field calculation is conducted to determine the phase dia¬ 
gram. For simplicity, we first consider the situation with 
M = 0 and (j) = tt/2 in in Eq. In our calcu¬ 
lation, we enlarge the unit cell to six sites of each hon¬ 
eycomb, and no further assumption for order parameter 
nii at these six sites are imposed, (i.e. totally 18 pa¬ 
rameters are determined from self-consistent iterations.) 
Enlarging the unit cell in the magnetic ordered phase 
turns out to be crucial for obtaining the state with lower 
energy and establishing the correct picture as discussed 
below. After we obtain the self-consistent solution, we 
can straightforwardly compute the single-particle gap for 
fermions, scalar spin chirality order and the Chern num¬ 
ber for mean-field ground state m- The resulting phase 
diagram is shown in Fig. 2, which contains both two sce¬ 
narios of phase transition, depending on the ratio ^ 2 /^ 1 , 
as well as four different phases: I. topological band insu¬ 
lator with no AF order (|m| = 0, C = 2); II. topological 
AF insulator with collinear AF order (|m| 7 ^ 0, C = 2 and 
5 = 0); III. trivial AF insulator with collinear AF order 
(|m| 7 ^ 0, C = 0 and S = 0), and IV. trivial AF insulator 
with canted AF order (|m| 7 ^ 0, C = 0 and S 7 ^ 0). 

First, for small 12/^11 such as the trajectory labelled by 
(a) in Fig. as t//ti increases, the system first undergoes 
a second-order phase transition across which a collinear 
AF order develops continuously (Fig. |^al)). As such 
a magnetic order increases, the mean-field Hamiltonian 
Hmf acquires a M {—M) term for spin-up (down), which 
suppresses the single particle gap at K [K') point (Fig. 
[^a2)). Thus, Hmf undergoes a trajectory as labeled by 
(a) in Fig. [^a). As |iTi| increases to a certain value, the 
single particle gap closes at K {K') point, beyond which 
the mean-field ground becomes a topological trivial one 
(i.e. C = 0, Fig. [^a3)). Along this trajectory, the spin 
chirality S is always zero (Fig. |^a4)). 

Secondly, for intermediate such as the trajectory 

labelled by (b) in Fig. as U/ti increases, the system 
first develops a collinear AF order (Fig. [|;bl)). Then, 
instead of reaching a topological transition, the system 
undergoes a first-order transition across which the mag¬ 
netic order becomes canted. This is accompanied by a 
jump of spin chirality order S (Fig. [^(b4)). This corre¬ 
sponds to a discontinuous change of effective M and5 in 
Hmf, as shown by trajectory labeled by (b) in Fig. j^a). 
Consequently, the system jumps from a topological phase 



FIG. 3: The magnetic order |m|, gap of fermion excitation 
A, Chern number C and scalar spin chirality order 5 as a 
function of U/t\ for three different value of as labeled 

by (a-c) in Fig. 

to a topological trivial phase (Fig. ib3)), and the gap 
closing point is avoided (Fig. [^b2)). This canted AF or¬ 
der can also be understood by a ring-exchange spin model 
in term of local spin picture. As U/ti further increases, 
the effective ring exchange is suppressed, S gradually de¬ 
creases and the system returns to a collinear AF insulator 
(Fig. §b4)). 

Finally, for large t 2 /H, such as the trajectory labelled 
by (c) in Fig. § as U/ti increases, a first-order transi¬ 
tion directly brings the system from a topological band 
insulator to a trivial canted AF insulator, across which 
|m|, C, S all display discontinuity (Fig. §cl),(c3),(c4)), 
and the fermion gap A remains finite all through (Fig. 
§c2)). 

Fluctuations Beyond Mean-field. Here we focus on 
the topological transition from Phase H to Phase HI, 
at which the gapless fermions coexist with gapless Gold- 
stone spin-wave mode of AF order. This invites the ques¬ 
tion whether the spin-wave fluctuation will change the 
critical behavior. To answer this question, we introduce 
the following low-energy theory with action 

S = J dtd^T (£„ +Ci + Cl) (5) 

^ [(i^tn)^ - c^(Vn)2] (6) 

£f = 4'^ [idt + ivFTzCTxdx -I- ivFdydy - toTzCTj,] vF (7) 

A =-A4''l'[cr2 G (n • s)]4'. (8) 

where £„ is a nonlinear sigma model that describes the 
low-energy fluctuation of AF Neel order n, and the eight- 
component object 41 = "^asa describes the Dirac fermion 
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FIG. 4: The spin-wave-fermion vertices are shown in (a,b), 
and the self-energy correction for the low energy fermions is 
shown in (c). The induced interaction among gapless fermions 
is shown in (d,e). The thin solid lines are the gapless fermion 
Cre and the thick solid lines are the gapped fermions Xk- The 
dashed line represents spin-wave ip = + iuy. 


nearby two valleys at K and K' points, in which s =t (i) 
denotes the spin up (down), a = 1,2 denotes the valley, 
and a = A{B) denotes the sublattice. The Pauli matrices 
Sx,y,z, '^x,y,z, and Tx^y^z refer to the spin, sublattice, and 
valley degrees of freedom, respectively. The parameters c 
and vf are the spin-wave and fermion velocities, respec¬ 
tively. Parameters g and A are coupling constants of spin 
fluctuation and coupling between spin and fermions, re¬ 
spectively. Here c, vp, g, m and A can all be given by 
microscopic parameters. In particular, in Eqj^the mass 
m is given by m = 3 -\/ 3 t 2 - 

With collinear AF order, we can assume that (n) is 
ordered along z direction, and then take riz = 1 and 
expand the action to the linear order of nx and ny, p = 
rix + ifiy representing a complex gapless boson field. To 
bring C into a more convenient form, we implement a 
local spin rotation 'I'(x, t) —)■ {7(x, t)4'(x, t), with U = 
exp[i{nySx — nxSy)l2 -!-••■], such that 17'^(n • s)U = Sz- 
Since this rotation involves only fermions, the Cn term is 
unchanged, while C{ + Ci becomes 

Cf + Cl — [iDt + ivpTzCTxDx + wpCTyDy 

-mazTz - XazSz] 4 *, ( 9 ) 

where the covariant derivative iDy^ = idy — ^{s+dyp* — 
S-dyp), in which sj- = ^{Sx^tiSy) and fj, = t,x,y. Thus, 
the mass term becomes —0 {rriTz (S> I + A/ 0 Sz)4', 
where I denotes 2x2 identity matrix. It is therefore 
clear that when X = m, namely, at the mean-field critical 
point, ^ki^ and Ik 2 t become gapless (the sublattice index 
a is suppressed). Hereafter we shall define = 'ki^ and 
^2 = 41 2 ^, which are gapless fermions, and xi = 'I'if and 
X 2 = 'k 2 ^, which are gapped fermions. 

The low-energy spin waves have small momentum, 
therefore the spin-wave-fermion interaction given by Dy 
terms does not change the valley index, i.e. by interact¬ 
ing with spin waves, can only turn into Xk with same 
K (Fig. s a,b)). After integrating out the spin-wave and 


the gapped fermions we can obtain a self-energy cor¬ 
rection for low-energy fermions as well as effective 
interactions among the massless fermions [12) . with cor¬ 
responding diagrams shown in Fig. |^d-e). To the lowest 
order the self-energy takes the form of Thus, it 

merely shifts the phase boundary. The induced interac¬ 
tion reads 


V = 
cf 


dPki d^k2 d?q 


(27r)3 (27r)3 (27r) 


3 ^. 


.k=1,2 


^ (^K,fc2-|-9^2;CK./C2) 2(^]^_j,^_^(T2:^iXi)(^2,fc2+9®'^^2,fc2)( lb) 


where q = (w, q), <Pq\s a, shorthand notation for J dqduj, 
and similarly for d^k. Neglecting the q dependence of 
V{q), we have V{q) = —u = — (1/2)'^g'^cA^ un¬ 
der certain approximation [T^ . where A is a momentum 
cutoff. 

This spin-wave-induced interaction, if sufhciently 
strong, can open up a gap at the nominal critical point. 
To see this fact, we only need to do a mean-field approx¬ 
imation of V. We find that if m > itc = ttvf/A, the gap¬ 
less “ground-state” at A = m is unstable towards the dy¬ 
namical generation [11 j of a mass term with 

A = Trvp{l/uc — 1/w). Away from the X = m point, the 
— sign (- 1 - sign) in ±A^^azSz^ is selected at A = to -I- 0+ 
(A = TO — 0+). That is to say, the generated mass jumps 
by 2A across the mean-field transition point A = to, thus, 
the nominal gap closing of fermions is avoided, and the 
transition becomes a first-order one. 

Finally, we remark that this physics triggered by mass¬ 
less spin-wave has no counterpart in the Kane-Mele- 
Hubbard model, because the SU{2) spin rotational sym¬ 
metry is explicitly broken there, thus the Goldstone mode 
is absent therein. 

Final remarks. Recent cold atom realization of the 
Haldane model can be naturally described by this HH 
model. In fact, in the experiment reported in Ref. [S], 
Mott insulator with suppressed double-occupancy sites 
has been observed. Our theoretical predications can be 
directly verified in this setup. In this realization, since 
the most crucial next nearest hopping term is generated 
by periodically shaking optical lattices, the periodic driv¬ 
ing will also modify the interaction term, in the order of 
l/oj (lo is shaking frequency). The thermal fluctuation 
of magnetic order may also be important. These effects 
will be left for future investigations. 

We would like to thank Yi-Zhuang You, Hong Yao, Fa 
Wang and Shou-Cheng Zhang for discussions. This work 
is supported by Tsinghua University Initiative Scientific 
Research Program, NSFC under Grant No. 1I304I75 
(ZW), No. 1I3254I8 (HZ), No. 11174176 (HZ), and 
NKBRSFG under Grant No. 2011CB921500 (HZ). 

Note added: Upon finishing this work, we became 
aware of Ref. m, in which the same model is studied. 
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SUPPLEMENTAL MATERIAL 
Mean-field theory 

The Hamiltonian of the Haldane-Hubbard model is given by 

H = -i2 E (e*^‘^4,A^ + h.c.) 

(b),s 

( 11 ) 

2,S i 

In the following we will only consider the M = 0 case. At the mean-held level, we decompose the on-site interaction 
term as 




ic77V - 

o q ^ 



-m.i ■ Sj -h 



( 12 ) 

(13) 
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where ^ s’^ss'Ci^s' is the spin operator and mi = AU {Si) /3 is the on-site magnetic order parameter. The 

ss' 

mean-field Hamiltonian can be constructed as 

^MF = -^1 H + h.c.^ - ta +h.c.) 



I 


This Hamiltonian is quadratic form and can be directly diagonalized. As explained in the main text, we enlarge the 
unit cell to a full hexagon containing six sites. We do a Fourier transformation 

= (R), 

R 

where R is the position of the unit cell, a = Al, Bl, A2, B2, A3, B3 denote the sublattices and Af is the total number 
of unit cells. The Hamiltonian can be transformed into momentum space: 


Hmf = d/t(k)Fr(k)d/(k), 


kGBZ 


where (k) is a 12-component spinor: 




st 




At 


At 


At 


At 


At 


At 


At 


At 


Diagonalizing the matrix H (k), one can obtain the energy band for the mean-field Hamiltonian 


Y,ulH,,{-k)U,, = S^,E^{k). 

u' 


(15) 


(16) 


(17) 


Therefore the diagonalized mean-field Hamiltonian is 

^MF= 5] (k)i?^(k)ci>^(k), 

kGBZ fJ, 


where (k) = ^ (k) is the fermion operator of each band. For the half-filling case, the ground state is the 

i 

full-filling of the lowest 6 bands. 


iGS) = n n (k) io) • 

^— 1 kGBZ 

The magnetization rria can be calculated from the ground state as 


m 


z 

OL 


m 


X 

OL 


m 


y 

OL 


^ (cft(’^)^“d(k)-c),_^(k)ca,;(k)^ 

kGBZ 

(cft(^)^^d(k)-cl.i(k)cc.t(k)) , 

kGBZ 

^ (eft Wcr.l(k)+c);, ;(k)c„,t(k))^g. 


(18) 


Since there are six sites in each unit cell, one has to self-consistently calculate 18 parameters of me,. We numerically 
calculate the ground state and iterate until the magnetization at each site converges. We then calculate the spin 
chirality order defined by 

5=(s.)-((s,)x(s,)), 

The Chern number of interacting quantum anomalous Hall insulators has been defined in Ref.[T]. At the mean-field 
level it is reduced to the Chern number of the mean-field wavefunction. We numerically calculate this Chern number 
using the algorithm given in Ref. [5]. 
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Effective Field Theory 


Near the phase boundary between Phase II to Phase III, we introduce the following low-energy thoery with action 

S = J dtdPr (£„ + £f + Li), 

^ {dtnf - (Vn)^ , 

Zg I J 

£f = [idt + ypT^a^id^ + vpayidy - 

Cl = —A'l'^ [az 0 (n • s)] 'b. 


The parameters are explained in the text. With the collinear AF order, we can assume that the n is ordered along z, 
direction, namely, ft « (0, 0,1). We expand £„ to the linear order of Ux and Uy as 


i—x^y 


(19) 


Then we make a spin rotation, 'I'(a;) —>■ [/(a:)'I'(x), with U{x) = exp [i (ny{x)sx — nx(x)sy) /2 such that 

U\x) (n • s) U{x) = Sz- 


Then the last term of the Lagrangian becomes 

Cl = {az ® Sz) ( 20 ) 

It is a mass term of the fermions, so it could be absorbed into the action of fermions. We have Cf + Ci 

Cl = [iDt -I- vpTzUxiDx + vpcryiDy - mTzCTz - XcTzSz] 'k, (21) 

where the covariant derivative is given by: 

Dfz = d^ + U^d^U 

= + ^ {s-dy,(p - s+df,(p*) , 

where s± = ^ (sx ± tsy), and (p = Ux + ifiy. This Ci + Ci can be written explicitly as 


Cl = [idt + vpCTxidx + ypUyidy - [m-f X) az] 'I'l-f + 'I'l^ [idt + vpcrxidx + ypayidy - (m - A) az] Ti; 

[idt - vpCTxidx + ypayidy + {m- X) az] 'I'2t + - ypaxidx + ypayidy -I- (m -b A) az] ^'2; 

+ ^^4 + ypc^xidxC + ypayidyip) (idtip* + vpaxidxip* + ypayidyip*) 

+ ^^4 - ypCTxidxC + vpayidyif) [idtC* - vpaxidxip* + ypayidyip*) T2;. (22) 

When A = TO, namely, at the mean field critical point, Ct becomes 

Cl = [idt + ypaxidx + vpayidy - 2maz] 'I'lt + 'I'l^ [idt + vpaxidx + vpayidy] 'I'l; 

-l-'I' 2 t [*'5* “ ypaxidx + ypayidy] ^ 2 ^ -b - ypaxidx + ypayidy + 2maz] 'I' 2 ; 

[idtip + Vpaxidxip + ypayidyip) 'I'l-^ - (idtip* + vpaxidxip* + ypayidyip*) 'I'l; 

+ ^^4®- - Vpaxidxip + ypayidyip) ^'2t - ^^45-1- {i^tip* - vpaxidxip* + vpayidyip*) ^<2;. ( 23 ) 

We can see that and 'I' 2 ^ become gapless, while and 'I' 2 ^ are gapped. We redefine 

Cl = 'I'U) C 2 = 'E' 2 ^, xi = X 2 = 'I' 2 ^ 

Cl = C 2 = 'i’l^crz, Xi = X 2 = 'i’licrz 

ip = nx + iriy, ip*=nx-iny, 



then the total action becomes 


= 

Ci = 

Cl = 


^ {\dM^ -C^ |Vv5|^) 


25 

6 (7o*9t + VFliidx + VFl2idy) 6 + 6 {loidt - VF^iidx + VFl2idy) 6 
+Xi [loidt + VFliida: + VFl 2 idy - 2m) xi + X2 - VFliid^ + VFl 2 idy + 2m) X2, 

iloidt^p + VFliidx^p + VFl2idy(p) Xi - ^Xi + VFliidx^* + VFl2idyV*) Ci 

+ ^X2 hoidt'f - VFliidxtp + VFX2idyip) 6 - ^6 {xoidt^* - VFXiidx(p* + VFl2idy(p*) X2- 

Note that we have mapped the low energy theory into four Dirac fields interacting with a complex scalar field. In the 
following, we will integrate out the massive Dirac fields and the complex scaler field to obtain an effective theory for 
the massless fermions. 

The partition function is given by 


Integrating out the x ip field gives rise to an effective action for ^ field 

giSeff ^ giS£ J (e*®0o ’ 

Here (e*'®^)p is the average over free x (f field, which reads 

= exp (Si) + — - 3 (Sf) ^ H ^ 

The effective action has the form of 


(24) 


(25) 


(26) 


(27) 


The second term in Eq. 27 generates a self-energy for the gapless fermion, which is illustrated in fig.4 (c). 


^ (^i^) = 


2 ! 


(p) (p) (p). 

( 27 r) 


This self-energy has the form of 


Si,2 (p) = z y 


d^k Xoko - VFllkl - VFX2k2 ■r^n\ 70^0 - VFXlkl - VFX2k2 . 


{2Fr 

■ d^k 


iD (k) 


(fcg - Vpkl - Vpkl) D (fc) iXi,2(p - k), 


-iKi,2(p- ^)> 


(28) 

(29) 


where D {k) is the propagator of the complex scaler field, and Ki^ 2 {k) is the propagator of the massive Dirac fermions: 

2g 


D{k) = 
-^1,2 (^) = 


1.2 _ ^21.2 _ ^2U2 I 
A/Q L- U rP2 

1 


le 


7ofco T VFXiki — VFl 2 k 2 T 2m -I- ie 
For low energy processes, we approximate K 12 (k) « so that the self energy is given by: 




d^k fcg — Vpki — v\k2 
(27r)^ fcg ~ ~ c^^2 + 


(30) 

(31) 

(32) 











9 


One can see that it merely shifts the mean-field phase boundary without qualitatively changing its physical properties. 

The third term in the effective action 27 as illustrated in fig.4 (d) and (e), generates an effective interaction between 
the gapless fermions: 


where the Vg is given by: 


1 f (Pq cPki cPk2 


Vg = 


i 1 


(27r)^ (27r)^ (27r)‘ 


d^k 


E ^n,ki —q^K,k2-\-q^K,k2^K,,ki ^^l,ki —q^2,k2+q^2,k2^1,ki 1 ; 


24 4^2 J (27r)^ 

X (/cq - Vpkl - Vpkl) D {k) {qo + k^f - v% {qi + ki) - Vp {q2 + k2) D {q + k). 


When g = 0, we have: 


- Vpkl) D'^ik). 


ig^ 

f 1 


2'^w? j 

(2^)^ 

\kQ — c^ki — c^/cl +ie ) 


[i _ /^i EV , A _ EV ± 

2‘^Tr^rrV 6 \ / 9 \ / 45 

where A is a momentum cutoff. This is the effective interaction between the fermions. 
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